One- and two-dimensional quantum walks in arrays of optical traps 
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We propose a novel implementation of discrete time quantum walks for a neutral atom in an array 
of optical microtraps or an optical lattice. We analyze a one-dimensional walk in position space, 
with the coin, the additional qubit degree of freedom that controls the displacement of the quantum 
walker, implemented as a spatially delocalized qubit, i.e., the coin is also encoded in position space. 
We analyze the dependence of the quantum walk on temperature and experimental imperfections 
as shaking in the trap positions. Finally, combining a spatially delocalized qubit and a hyperfine 
qubit, we also give a scheme to realize a quantum walk on a two-dimensional square lattice with the 
possibility of implementing different coin operators. 

PACS numbers: 03.67.-a, 32.80.Pj, 42.50.-p 
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I. INTRODUCTION 

In classical computation, random walks are powerful 
tools to address a large number of problems in many ar- 
eas of science, as, for example, graph-connectivity or sat- 
isfiability problems 0] . It is this success of random walks 
that motivated to study their quantum anàlogues in order 
to explore whether they might extend the set of quantum 
algorithms. Two distinct typcs of quantum walks have 
been identified: for the continuous time quantum walk a 
time- independent Hamiltonian governs a continuous evo- 
lution of single particle in a Hilbert space spanned by the 
vertices of a graph while the discrete time quantum 
walk rcquires a quantum coin as an additional degree of 
freedom in order to allow for a discrete time unitary evo- 
lution in the space of the nodes of a graph. The connec- 
tion between both types of quantum walks is not clear up 
to now 3j , but in both cases different topologies of the 
underlying graph have been studied, e.g., discrete time 
quantum walks on circles 0, on an infinite line p|, on 
more-dimensional regular grids and on hypercubes 
Q. The field has recently been reviewed by Kempe Q. 

Several algorithms based on quantum walks have been 
proposed H,H, El El . To implement such an algorithm 
in a physical system it ultimately has to be broken down 
into a series of gates acting on a register of qubits Q. 
From the more fundamental point of view however more 
straight-forward implementations are interesting, i.e., di- 
rect implementations of a quantum walker (a particle, 
a photon etc.) moving, e.g., in position or momentum 
space. So far some setups for one-dimensional realiza- 
tions have been analyzed, including trappcd ions |l2j| . 
neutral àtoms in optical lattices with state dependent 
potentials [ï^, single-photon sources together with lin- 
ear optical elements |Ï4|, and also with classical òptics 
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|l5| . Here we use the idea of spatially delocalized qubits 
(SDQ) developed in |16| to propose a novel quantum walk 
implementation with neutral àtoms. The particle is walk- 
ing in position space, but in contrast to the proposal in 
|l3| also the quantum coin is represented by a spatial de- 
gree of freedom, as it is implemented by the presence of 
the atom in the ground state of one out of two trapping 
potentials. The particle is manipulated only by varying 
the trapping potentials, which induecs tunneling between 
traps, and no state dependent potentials are necessary. 
This concept can be applied to neutral àtoms trapped 
in optical lattices 0, in magnètic potentials 0, as 
well as in arrays of microtraps |Ï9j : here we will espe- 
cially analyze the latter case. We will also show how a 
combination of a spatially delocalized qubit and a hyper- 
fine qubit together with state dependent potentials al- 
lows to implement a quantum walk on a two-dimensional 
square lattice. Quantum walks in higher dimensions offer 
a very rich structure of dynamics, and recently a spatial 
search algorithm using a modified quantum walk on a 
two-dimensional grid has been proposed [ïïj . 

For the one-dimensional case we will discuss the in- 
fluence of non-adiabatic processes and of shaking of the 
trap positions, and we will estimate the cffcct of decoher- 
ence. We will also consider dependències of the quantum 
walk on the vibrational trapping state and thus on the 
temperature and show that, within a range of parameters 
accessible in experiments, a transition from the quantum 
walk to the classical random walk can be studied. This 
is not only interesting from a fundamental point but also 
allows to assess the degree of control that can be reached 
in the experiment. 

It has been noted 0, |2(], Ü0|j that cssentially only 
interference is necessary for a quantum walk, such that 
it can be implemented with classical fields. Neverthe- 
less considering setups with neutral àtoms is justified by 
a strong interest in these systems as tools for quantum 
computation [2^ , as well as by the possibility to include 
further cffects as, e.g., quantum walks with two or more 
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(possibly interacting) partides [23l| . 

II. QUANTUM WALKS AND OPTICAL 
MICROTRAPS 

1. Quantum walks 

Let a particle move on a one-dimcnsional infinite linc, 
such that it can only hop between discrete sites x = k ■ a 
labeled by k £ {..., —2, —1, 0, 1,2,.. .}, with a being the 
distance between sites. At each time step the parti- 
cle moves with equal probability to either of the adja- 
cent sites. For a classical random walk, the probabil- 
ity for the particle to be at a certain site for a large 
number of steps approaches a gaussian function cen- 
tered around its initial position xo, with the variance 
(T 2 = ((x — .to) 2 ) growing linearly with the number of 
steps n. For the quantum version a state |fc) is attached 
to each site x = k ■ a, i.e., the particle is walking in 
Hw = span{|fe), k = . . . , — 2, — 1, 0, 1, 2, . . .}. However, 
the random move cannot be just replaced by walking to 
the left and to the right in superposition, as this turns 
out to be non-unitary |24j . For this reason a quantum 
coin is introduced as an additional degree of freedom. 
In the simplcst casc of the quantum walk on a line, the 
coin space Tic is two-dimensional and we will denote the 
states that span Tic by |— ) and |+), and the total Hilbcrt 
space is thus TLw ® 7~Lc- Each step of the quantum walk 
is then composed from two operations: (i) applying a 
unitary operation C to the coin (simultaneously at all 
sites), e.g., a Hadamard operation C = H: 

(I®H)|*,±) = -^(|fc,+>±|fc ) ->) Vfc, (1) 

followed by (ii) applying a displaccmcnt operation O^d 
which moves the particle left or right depending on the 
coin: 

Oid|*,±) = |À!±1,±) Vfc, (2) 

wherc we have not cxplicitly written the tensor product: 
\k) ® |±) = |fc, ±) ctc. The probability distribution aris- 
ing from the iterated application of W = Oid(I ® H) is, 
except for the first three steps, significantly different from 
the distribution of the classical walk: if the coin initially 
is in a suitable superposition of |— ) and |+) it has two 
màxima symmetrically displaccd from the starting point. 
In general the exact form of the distribution, especially 
the relative height of the màxima, depends on the ini- 
tial coin state. Compared to the classical random walk 
its quantum version propagates faster along the line: its 
variance grows quadratically with the number of steps n, 
a 2 oc n 2 , compared to a 2 oc n for the classical random 
walk. 

For the walk on a line, H is, up to phascs (which can 
be absorbed also into the initial state), the only unbi- 
ased coin operator pEj . For a two-dimensional regular 



square lattice a much richer structure of coin operators 
and possible probability distributions arises. As has been 
observed by Mackay et al. j|| and by Tregenna et al. [2í| 
in this case different unbiased coin operators and initial 
states can be chosen that produce significantly different 
dynamics, ranging from distributions with a sharp cen- 
tered spike to distributions having the shape of a ring. 

2. Optical microtraps 

As a particular setup for the implcmcntation, we con- 
sider the controllcd motion of neutral àtoms in arrays 
of optical microtraps. The microtraps are created by il- 
luminating a set of microlenses with a red detuned làser 
beam, such that in each of the foci of the individual lenses 
neutral àtoms can bc stored by the dipole force 0]. 
illuminating the set of microlenses by two independent 
làser beams, it is possible to generate two sets of traps 
which can be approached or separated by changing the 
angle between the two làsers. This allows the atom to 
propagate between different microtraps. The optical po- 
tcntials have a gaussian shape, i.e., 

V{x) = -Vbcxp |-^™iï 2 | = -Vbexp (^--j^r( ax ) 2 

(3) 

For the simulations we present here, we will use Vq = 
200 huj x . In this case the traps are deep enough to be 
described by harmònic potcntials of frequeney lo x in the 
limit of large separation. Then a^ 1 = yl h/muj x denotes 
the spread of the ground state in position space, with m 
being the mass of the atom. 

For the preparation of the initial state we assume that 
a single atom can be placed in the ground state of a 
specific trap. For this reason, and also to be able to read 
out the final state of the system, it is necessary to be 
able to address each trap separately. This addressability 
has already been demonstrated in arrays of up to 10 x 10 
microtraps |l9j . 

III. ONE-DIMENSIONAL WALKS 

The implementation of the coin at each site k will fol- 
low the idea of spatially delocalized qubits from i.e., 
the basis states |fc, ±) will be represented by a single atom 
occupying the ground state of one of two adjacent traps. 
Unitary operations are performed by approaching the two 
traps forming the coin, allowing the atom to tunnel be- 
tween them. In the following we will use quantum òptics 
notation to describe the effect of tunneling between traps, 
e.g., an operation exchanging the population of two traps 
will be termed tt pulsc and a Hadamard-likc operation 
|fc,±) i— > ± i\k, — )) will be termed ir/2 pulse 

We propose two closely related configurations, both 
lcading to a quantum walk. For the first configuration 
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FIG. 1: (Color online) Configuration with two rows of traps, 
the qubit is implemented 'perpendicular' to the rows (dashed 
rectangles show which two traps form each qubit). The upper 
(lower) row moves left (right) with constant velocity (see ar- 
rows on right). (I) After the firat step: \tp) = l/y/2(\ - 1, -) + 
| + 1,+)); (II)- (III) thecoin operation, inthis case a Hadamard 
gate, is performed when the traps pass each other; (IV) the 
shift Om is implicit through a redefinition of the qubits. After 
an even (odd) number of shift operations only the even (odd) 
qubits (compared to the Standard quantum walk definition) 
are defined; (V) the probability distribution after the sixth 
displacement operation. For the numerical simulation we used 
a potential which, along the line connecting the centers of 
two traps, reads V(x) = hw x min {a 2 (x — a) 2 , a 2 (x + a) 2 } . 
The velocity is chosen such that during the passing of two 
traps a Hadamard operation is performed. The initial state 
is |V'init) = -±5(|0,-) + |0,+». 
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one row with respect to each othcr rcquircs mcchanical 
movement of an array of lenses this setup is quite chal- 
lcnging, a problem which might be overcome by using 
holographic tcchniques to generate arrays of microtraps 

m 

In what follows we will concentrate on the coin being 
implemented 'paral·lel' to the direction of displacement, 
such that only a single line of traps is necessary, see Fig-El 
Labeling the traps of the fcth qubit by 2k and 2k + 1, for 
coin operations the traps 2k and 2k + 1 are approached, 
whilc for the steps in the walk a 7r-pulse between traps 
in adjacent qubits, i.e., between traps 2k + 1 and 2(k + 
1), moves the atom one step to the left or to the right, 
respcctively. Contrary to the displacement opcrator from 
Eq. J5J), this procedure flips the coin operator at each 
move, i.e., we have Om\k,±) = \k ± 1,+) (termed flip- 
flop walk in ^3). Clearly the experimental requirement 
is to be able to move all odd (or all even) traps as a 
whole to both directions, thus approaching each second 
trap to its left or right neighbor. This can be realized as 
described in section III 2 1 in optical microtraps ; but 
also in optical lattices |27l l28j or magnètic microtraps 
0. 
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two rows of traps arc necessary. Each coin is defined 
through one trap from each row. By moving both rows 
in opposite directions with appropriately chosen distanec 
and velocity, the coin operations are performed when the 
traps pass each other at close distance. The displace- 
ment is implicit through a redefinition of the coin each 
time two traps have passed. Fig.^ shows the first steps 
in the temporal evolution of the configuration with two 
rows of traps along with the corresponding probability 
distributions resulting from an integration of the two- 
dimensional Schròdingcr cquation (see figuro caption for 
details). Fig. Q] (I-III) show the coin operation, (III- IV) 
the redefinition of the coins and (V) shows the probability 
distribution after the sixth displacement operation. The 
onset of the quantum walk character of the distribution 
is clearly visible as two màxima symmetrically displaccd 
from the origin appear. If the continuous displacement of 



FIG. 2: (Color online) Configuration with one row of traps, 
the qubit is implemented 'parallel' to the rows (grey boxes). 

(I) After the first step: = (| - 1,+) + | + 1, -»/>/2; 

(II) traps inside each qubit are approached to give the coin 
operation; (IV) the shift Oid is realized through approaching 
traps of adjacent qubits. 

For the gaussian trapping potentials of Eq. (J3J) , Fig. [3| 
(i) shows a simulation of the quantum walk for po- 
tential depth Vo = 200 fko x and an initial separation 
Qíamax = 60, obtained from an integration of the one- 
dimensional Schròdingcr equation using Fouricr transfor- 
mation and a split-stcp method. Initially the atom is pre- 
pared in an equal superposition of the two ground states 
of the central qubit, i.e., of the two central traps, such 
that |Vw) = -75(10,+) + |0,-)). The distanec of the 
traps is changed between the maximal value aa max = 60 
and a minimal value aa mm = 28.8. The latter distance 
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is for the givcn trapping paramcters close enough for 
tunneling to take place. Moving the traps adiabatically 
between this distances requires techniques to optimizc 
the moving process while suppressing transitions between 
motional states |29| . In this way, the time t r necessary 
to approach - or separate - the traps can be reduced to 
Lu x t r = 100 or below while maintaining a fidclity largcr 
than F = 0.999 \Ü^. The time U for which the traps 
are kept at the distance a m i n is chosen such that alter- 
nately a 7r pulse and a ir/2 pulse are applied. The figure 
shows the population of the traps after t = 10, 20, and 
25 steps. In Fig. [31 (i) the characteristic shapes of the 
quantum walk distributions are visible. Subsequently we 
will analyse how the probability distribution changes if 
differcnt vibrational states are involved or experimental 
impcrfcctions are present. 

1. Excited vibrational states: the influence of temperature 

Tunneling as well as adiabaticity do crucially dcpend 
on the timing of the change of the trap separation. For 
all simulations t r , the time needcd to move the traps 
together or apart, and Í£, the time for which the trap 
separation is kept constant, are chosen to apply the cor- 
rect operations for the vibrational ground state. If the 
atom starts in an excited vibrational state, then the tun- 
neling rate is larger and thus in general the coin operator 
C as well as the displacement operator Oid change. The 
former will be distinct from the Hadamard operator H 
and in general biased, 

C'\k,±) = y /p\k,+)±jï=p~e iA °\k,-), (4) 

(the Standard unbiased Hadamard operator has y = \ 
and Ac = 0), the latter will take a general form 

0' 1D |fc,±) = V^|fc±l,T>±yi _ ^e ,;Ao |fc,±). (5) 

(c = 1 and Ao = for the Standard displacement oper- 
ator). For an atom in a fixed vibrational level, p, Ac, c 
and Ao, and thus the operators C and 0' 1D are con- 
stant, because the movement of traps is assumed to bc 
unchanged throughout the process. In such a case the 
qualitative shape of the probability distribution is not 
modificd significantly, it still shows the characteristic 
symmctrically displaced peaks. However, a simulation 
for an atom initially in the first excited vibrational state, 
c.f., Fig.0(ii), shows a distribution which csscntially has 
a central peak and long symmetric tails. In this case 
the variance a 2 grows only lincarly with the number of 
steps, as compared to a quadratic increase of the vari- 
ance for the ground state distribution. The difference 
to the expected result can be attributed to the fact that 
the approaching and separating processes were optimized 
to suppress non-adiabatic excitations from the ground 
state. For highcr vibrational states excitations are non- 
negligible, causing coin as well as displacement operator 
to induce transitions between different trapping states. 



Then effcctively wc have a one-dimensional walk with 
a highcr dimensional coin. A quantum walk distribu- 
tion should be reobtained when restricting the quantum 
walk to some fixed vibrational state by suppressing non- 
adiabatic transitions. This can be done by increasing the 
time t r used to approach the traps. Then, as can be seen 
in Fig. |3 (iii), again the characteristic displaced peaks of 
the quantum walk probability distribution appear. As 
for the ground state, the variance increases quadratically 
with the number of steps. Note however that the dis- 
tribution is not the same as for the ground state due to 
different coin and displacement operators. 

A more realistic assumption than starting from a pure 
state with the atom being in a specific vibrational level 
is to consider a thermal Boltzmann distribution of the 
vibrational modes 

oo oo 

p= lj2e-^\m, z^e-f**, fi=j±=, (6) 

3=0 3 

where Ej is the energy of the jth vibrational mode. In 
this case the experimentally accessible probability distri- 
butions are the classically averaged probability distribu- 
tions, weighted with factors exp(—(3Ej)/z. Rcspective 
probability distributions after n = 20 steps are shown in 
Fig. ^] for initial ground state populations of 50% and 
25%, corresponding to a mcan number of vibrational 
quanta of (i>) = 1 and (y) = 3, or to a temperature 
of T = 1.1 fjK and T = 2.7 pK (for Rb àtoms and trap 
frequeney u> x = 10 5 s _1 ), respectively. The character- 
istics of the quantum walk remain visible even at such 
(u). In optical lattices with parameters similar to what 
we consider here, ground state populations of above 98% 
have been achieved [3(| • Thus we can expect that the 
range of temperatures necessary to observe the quantum 
distribution is well within the reach of experiments. 



2. Experimental imperfections and decoherence 

As an important experimental imperfection we will an- 
alyze shaking of the centers of the traps. We assume that 
the movement of traps in the even or odd sets (traps 2k or 
2fc + l, respectively) is correlated. For an implcmcntation 
with optical microtraps this is justified, because each set 
can be generated from a single làser beam, as described 
in section III 21 but it will also describe the situation 
for optical lattices with tunneling controllcd by changing 
the intensity of one out of two counterpropagating làser 
beams, as in [27], [2^]. We also anticipate shaking with 
a fixed frequeney far away from the trapping frequeney. 
This will not eliminate transitions between vibrational 
states, as it alters the (optimized) path to approach the 
traps. Due to the strong sensitivity of tunneling on the 
distance, shaking will give rise to changes in the rate of 
the population oscillation between the traps, i.e., coin 
and displacement operators will change. They will even 
bc different from step to step, as shaking is not correlated 
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FIG. 3: The probability distribution to find the atom at a specific trap site, with (I) the ground state and (II) the first 
excited state as the initial vibrational state, for a one-dimensional quantum walk on a finite line of 62 traps; from top to 
bottom distributions after t = 10, 20, and 30 steps are shown. Parameters: Vó — 200fàcj^, aí<x max = 60, ctíimm — 28.8, and 
iL> x t r — 100; w x ti,Tv = 20.25 for the n pulse and w x í ij7r / 2 = 112 for the n/2 pulse (for simplicity we fixed a m i n and then searched 



for the smallest ti that produces the desired operation. As tunneling already happens for a > a n 



t, 



does not give the 



identity operation, and for this reason í; j7r /2 > íi.vr). 
are suppressed. 



(III) like (II), but with u x t r = 200, such that non-adiabatic excitations 




FIG. 4: Probability distributions after t = 20 steps for 
thermal Boltzmann distributions of vibrational modes, ini- 
tial ground state population (I) 50% and (II) 25%. All other 
parameters as in Fig. [3] 



to the global motion of approaching the traps. The con- 
sequences of this random variation should be similar to 
the effects observed in the presence of decoherence. De- 
coherence in quantum walks has becn studied by Kcndon 
and Tregenna in a general framework |3lj| and by Dür et 
al. for the spccial case of an optical lattice implcmcnta- 
tion of a quantum walk |l3j . In the presence of decoher- 
ence, the probability distribution of the quantum walk 
on the (infinite) line ultimately collapses to the gaussian 
distribution characteristic for the classical random walk. 
However, for the product of the numbcr of steps t and 
the decoherence rate p being small enough, decoherence 
does not significantly degrade the quadratic sprcading of 



the walk. As has furthermore been found in |3Ï|. for a 
certain intermediate choice of tp a highly uniform distri- 
bution of the probabilities between positions ±t/y/2 can 
be observed if decoherence acts on the position degree of 
freedom or on both, the coin and the position degree of 
freedom; note that for our implcmentation the effect of 
shaking corresponds to the latter case. 

In Fig. the results for a sinusoidal variation of the 
trap distances around the perfect value with frcqucncy 
^Shake = 0.01 üj x and amplitude aAa are shown. The 
transition from the quantum to a classical distribution 
takes place for amplitudes on the order of a percent of 
the minimal distancc, the intermediate flat distribution 
is clearly visible at aAa m 0.09. For larger amplitudes 
of shaking the non-adiabatic transitions are dominant 
(Fig. (II)), and the variance decreases strongly with 
increasing shaking (Fig. (III)). For smaller amplitudes 
however the variance initially increases with increasing 
amplitude of shaking. To quantify the flatness of the dis- 
tribution we also calculate the total variational distance 
v(t) = ^2 n \P(n,t) — P u (t)\ to the uniform distribution 
P u (t) of half width t/y/2 0. Here P{n,t) is the prob- 
ability to find the particle in the traps belonging to the 
n-th qubit after t steps. As Fig. [ïï] (IV) shows, the total 
variational distance decreases initially, before it increases 
again as the probability distribution approaches a gaus- 
sian. 

Let us discuss the duration of the operations neces- 
sary for the quantum walk in order to estimate the in- 
fluenec of other decoherence mechanisms in the experi- 
ment. As the processes rely on tunneling, the duration 
of a singlo operation is on the order of the inverse trap- 
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FIG. 5: The effect of shaking of the trap position on the 
quantum walk. Shaking is modeled by a sinusoidal variation 
of the trap distances around the perfect value, with frequency 
uishake = 0.01u x and amplitude aAa (all the other param- 
eters are as in Fig. [3J . (I) The probability distributions for 
various vàlues of aAa (after tracing over the coin degree of 
freedom). (II) Ground state population, (III) variance, and 
(IV) variational distance v from the uniform distribution af- 
ter í = 17 steps. In (II) and (III) dashed lines and squares 
give the respective vàlues for the full population, dotted lines 
and circles for the ground state only. 



tudc fastcr, as the initial separation of the àtoms can be 
made shorter, such that decoherence rates similar to the 
case of optical microtraps can be expected. In both cases 
these decoherence mechanisms affect the coin as well as 
the position, and for this case the crossover from the 
quantum to the classical distribution has been numeri- 
cally estimated in j3j| to take place at tp m 2.6. For this 
reason it should be possible to observe the quantum walk 
in such systems and to analyze changes caused by tem- 
perature and shaking without being limited by decoher- 
ence from photon scattering etc. The strong dependenec 
on temperature and on non-adiabatic transitions of the 
quantum walk with delocalized qubits might thus be in- 
teresting as a tool to analyze to which extent the ground 
state population, the shaping of the trapping potcntials, 
and tunneling processes can be controllcd for a particular 
experimental setup. In addition, quantum walks in this 
particular physical system could be used to investigate 
how decoherence acts with respect to the spatial degree 
of freedom. 



IV. A TWO-DIMENSIONAL QUANTUM WALK 

For the two-dimensional quantum walk on a regular 
square lattice, i.e., if 

H w = sp&n {\(k,l)),k and l e {...,-2,-1,0,1,2,...}}, 

a four dimensional coin degree of freedom to control the 
displacement of the particle into the four possible direc- 
tions is necessary: 



Tic = span{| + +) = 



/>}• 



(7) 



Here we propose to implement such a coin by a suitable 
combination of a spatially delocalized (SD) qubit and 
a hyperfine (HF) qubit combined with spin-dependent 
transport jïj, y3j, i-e., Tic is a tensor product of the 
Hilbert space formed from the ground states of two ad- 
jacent traps and from two hyperfine states of the atom: 
Tic = TÍsb ® Whf- 



ping fr equency, which typically is about ui x = 10 5 s _1 
|19l l32| . For the parameters used here a single applica- 
tion of Oro (I <8> C) takes around 5 ms. The dominant 
decoherence mechanism can be expected to be the scat- 
tering of photons from the trapping làser, with scattering 
rates on the order of 0.1 . . . 1 s" 1 0,H2- Then the prob- 
ability for a decoherence event to oceur within a single 
step is p = 0.0005 . . . 0.005 and for t = 17 applications of 
Om (I ® C) we have tp = 0.0085 . . .0.085. For optical 
lattices the same decoherence mechanism is present, but 
also decoherence through fluctuations in the phase of the 
làsers producing the lattice, giving rise to fluctuations 
in the trapping potentials should be taken into account. 
One the other hand operations can be an order of magni- 



1. Separable walk 

There is no unique extension of the Hadamard oper- 
ator H to Tic from Eq. Q, becausc diffcrcnt classes of 
unbiased coin operators for two-dimensional walks exist 
|25j . The most obvious and simple generalization is to 
take a Hadamard coin for both directions. This can be 
realized by first approaching the traps to perform a 7r/2 
pulse for the delocalized qubit as above, and then putting 
the atom in a superposition of the two hyperfine levels 
by a 7r/2 two-photon |34| or microwave 01 pulse, which 
realizes C2D = Hhf^Hsc ( c -f- Fig. El(a-c) for the case of 
| + +) as initial state). For the coin-dependent displace- 
ment assume that at each vèrtex of the two-dimensional 
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FIG. 6: (Color online) Implementation of the coin and the coin operator C20 = Hsd C3> Hhf for the two-dimensional walk: (a) 
The four levels are formed as a tensor product of a delocalized qubit (left and right traps) and a hyperfine qubit (dark and 
grey filled circles symbolizc the +} and — ) hyperfine states, full and dashed lines denote the respective trapping potentials). 
The initial statc is | + +). (a)— >(b) For Hsd the traps are approached (for both hyperfine states). (b)— »(c) For Hhf a n/2 
pulse between the two hyperfine levels is applied to all traps simultaneously. (d) The implementation of the two-dimensional 
walking operator O20 as a combination of tunneling and spin dependent transport: (A) traps in horizontally adjacent traps 
are approached to give a 7r-pulse as in the one-dimensional walk; (B) the lattice is displaced in opposite vertical directions for 
the two hyperfine states. 



grid each two traps forming a coin arc aligned horizon- 
tally. Then in horizontal direction first the walking op- 
erator Oid can be applied, i.e., within each row traps of 
ncighboring qubits are approached as described above to 
give a 7r pulse, followed by transiating the lattice poten- 
tial in opposite vertical directions for each spin state, as 
proposed in 0| ( see Fig. El (d))- In total, the action of 
the walking operator 2 d m 7~(-w ® 7~Lc is given by 



2D |(M),±±) = \(k±l,l±l),T±) 



(8) 



Fig.Hshows a probability distributions arising from al- 
tcrnatingly applying I ® C 2 d and 2 d to the initial state 
iVWt) = 1(0,0), ++). From its construction it is easy to 
see that the coin operator C 2 d does not mix the horizon- 
tal and vertical direction. For this reason one recovers 
the one-dimensional quantum walk when projecting the 
distributions along the x or y directions (due to the choicc 
of the initial conditions the distribution is not symmetric 
in this case). C 2 d is thus a separable Hadamard walk 
according to the classification of @ . 



2. Entangled walks 

Morc sophisticated coin operators are also possible, 
and we will show how to implemcnt one which entan- 
gles the two directions. Assuming to be able to change 
the trapping potentials for both hyperfine states indepen- 
dently, we can apply a 7r/2 pulse (Hadamard operation) 
on delocalized qubit for the |+) hyperfine state and a ir/2 
pulse followed by a tt pulse (Xnot operation) on the de- 
localized qubit for the |— ) hyperfine statc. Subscqucntly 
a 7r/2 pulse is applied tothe hyperfine qubit as beforc. 



Then, defining the Xnot operation as, 
Xnot = 



1 

1 



(9) 



the full coin operator reads 



r< 

L - 



Ent 
2D 



(I 


3H HF )-( 


, ySD ttSD 
"■" NOT 




/l 1 1 


1 


1 -1 1 


2 


1 1 - 
\l -1 - 




(10a) 
(10b) 



(Similary, given local addrcssability, diffcrent làser pulscs 
could be applied to the two traps forming the spa- 
tially delocalized qubit). The operator Cf^ is non- 
separable @,[2í|. The result for the initial state IV'init) = 
|(0, 0), ++) is shown in Fig. |SJ Similar sequences of op- 
erations can be used to g enerate different coin operators. 
It has been noted in [25j that even for a fixed entangling 
coin operator very different probability distributions can 
be obtained through varying the initial state. At least 
for a system of optical microtraps it should be possible 
to engineer the initial state carefully enough, because due 
to the large separation of traps, single sites can easily be 
addressed. In this way these systems can be an interest- 
ing testbcd to cxplore the rich structurc of probability 
distributions of two-dimensional quantum walks. 

Recently Ambainis et al. proposed a quantum search 
algorithm to locate a single marked site among N loca- 
tions arranged on a square V^V x s/N grid using a to- 
tal number of 0(VN~ \og 2 N) steps, thus outperforming 
Grovers algorithm in this case [ll|, where the time to 
move between locations is taken into account. The algo- 
rithm is bascd on a quantum walk on a two dimensional 



8 



(I) 



PÍU) 




Y trap index 1 



trap index k 



(II) 




j trap index I 



lattice with (i) periòdic boundary conditions, (ii) a spe- 
cial coin as well as a special displacement operator, and 
(iii) a certain initial state. We will show how to realize 
the bàsic ingredients (ii) and (iii) within our proposal. 
The realization of periòdic boundary conditions however 
is not straightforward for the setup proposed here, but 
might be achievable with more advanced setups. Fixed 
boundaries can be imposcd for optical microtraps by only 
illuminating a rectangular subset of lenses and apply- 
ing an additional ir/2 pulse on the hyperfine qubits on 
the borders of the lattice. Numerically we observe that 
the effect of a large amplitude at the marked site is still 
present for such boundary conditions, though the effect 
on the performance of the algorithm necds further invcs- 
tigation. 

The required displacement operator is the two- 
dimensional version of the flip-flop walking operator 
which we already saw in the one-dimensional realization. 
It reads 

Off|(M).±±) = \(k±l,l±l),TT), (11) 

i.e., the walk changes direction after each step, and obvi- 
ously here it can be realized by a 7r-pulse on the hyperfine 
qubit after applying C^d- The coin operator has to be 
chosen as 



Cn = 




(12) 



(HI) 



0.04 



P(k,l) 




trap índex k 



except for the marked vèrtex, for which 



d = 



-1 














-1 














-1 

















1/ 



(13) 



is applicd. For an implcmentation of the scarch algo- 
rithm, these operators have to be constructed as a suit- 
able combination of operations on the hyperfine and the 
delocalized qubit. We will require the manipulation of 
the delocalized qubit to act identically on all sites, but in 
order to engineer a special coin operator for the marked 
vèrtex, for the manipulation of the hyperfine state we 
will assumc to be able to address single sites. We definc 
a single qubit phasc gate 



-^PHASE 



1 

-1 



(14) 



FIG. 7: (Color online) Probability distributions for the two- 
dimensional quantum walk obtained from alternatingly ap- 
plying C*2d and O2D to the initial state IV'init} = 1(0,0),++) 
after (I) t = 1, (II) t = 2, and (III) í = 25 steps. The co- 
ordinates on k and l axes label the traps, thus in k direction 
each two traps form one site; note that the total population 
of each trap is shown (sum of probabilities for both hyperfine 
states) . 



and 



^NOT 



X 



HF 

NOT 



I, 



(15) 



as an operator which produces a tt pulse only on one trap 
of each coin. Thcn 



Cn — X[ 



NOT ' 



H 



SD 



X 



HF 

PHASE 



NOT ' 



(H 



SD 



^PHASe) ' -^NOT 



(16) 
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produces, except for an overall phase, the correct coin 
operator. The operator for the marked vèrtex is obtaincd 
by merely replacing the A^qt operator by I: 

-C x = (H SD (g) ^phase) ' (H SD <8> ^phase) • (l'O 

Thus, given the possibility to locally manipulate the hy- 
perfine qubit and except for a global phase, the quan- 
tum search coin operators can be constructed. Finally, 
the initial state has to be chosen as an eigenstate of 
Off • (I® Co). For fixed boundaries, such an eigenstate 
is given by 

i $ o> = 77^EE E E k*,o.«j9). 

ZVJV k=x i=i ae{+,-} 

which can be generated by a sequence of shift operations 
firstly via tunneling and secondly via displacements of 
the state-dependent lattices [ïïj . 

V. SUMMARY 

We have discussed the implementation of quantum 
walks with a neutral atom trapped in the ground state of 
optical potentials by using the concept of spatially delo- 
calized qubits, i.e., a coin dcfhied through the presence of 
the atom in one out of two trapping potentials. We have 
shown that in this case a quantum walk on a line can 
be performed in a simple way only through a variation 
of the trapping potentials, without the need for addi- 
tional làsers to addrcss intcrnal states of the atom. Our 
simulations were performed with realistic parameters for 
present optical microtrap systems, but the concept is as 
well applicablc to optical lattices or to magnètic micro- 
traps. We have studied the influences of various experi- 
mental impcrfcctions on the probability distribution and 
have found a strong change if the atom is initially not 
in the ground state of the trap. This change, leading 
to a strong dependence of the quantum walk on tem- 
perature, can be attributed to non-adiabatic excitations 
to other vibrational states during the movement of the 
traps. We have also studied the influcncc of shaking and 
found a transition from quantum to classical probability 



distributions, taking place for shaking amplitudes on the 
order of 1% of the tunneling distance. As an intermediate 
step, this transition exhibits a very flat distribution. An 
cstimate of other decoherence effects such as scattering of 
photons from the trapping làsers suggests that quantum 
walks should be observable in the experiment and the 
effects of temperature and shaking should be accessible 
to experimental investigation. In this way, implcmcnting 
the quantum walk with spatially dclocalizcd coins could 
give information on the extend to which the ground state 
population and the movement of the traps can be con- 
trolled. 

Finally, we have combined the concept of the spatially 
delocalized qubit with a hyperfine qubit and state depen- 
dent potentials to obtain a scheme to implement a quan- 
tum walk on a two-dimensional regular latticc. Within 
this scheme, which again is close to what is realizable 
with state-of-the-art technology in optical microtraps as 
well as in optical lattices, different coin operators are 
possible, such that in this setup the variety of different 
distributions in two-dimensional quantum walks can be 
explored. Especially we have shown how to construct 
separable and entangling coin operators, as well as the 
operators necessary to implement a spatial search algo- 
rithm on a two-dimensional grid. It is worth stressing 
that the scheme proposcd by us can be used to construct 
the generalized coined quantum walk, in which the walker 
aequires at each step a phase [3^ • 
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p(k,l) 




FIG. 8: (Color online) Probability distribution of the two- 
dimensional quantum walk with the coin operator Cfo* after 
t = 25 steps. The initial state is ji/'mit} = |(0, 0), ++). 



